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Abstract Chiral anomalies have profound impact on the transport properties of rel- 
ativistic fluids. In four dimensions there are different types of anomalies, pure gauge 
and mixed gauge-gravitational anomalies. They give rise to two new non-dissipative 
transport coefficients, the chiral magnetic conductivity and the chiral vortical con- 
ductivity. They can be calculated from the microscopic degrees of freedom with the 
help of Kubo formulae. We review the calculation of the anomalous transport co- 
efficients via Kubo formulae with a particular emphasis on the contribution of the 
mixed gauge-gravitational anomaly. 



1 Introduction 

Anomalies in relativistic field theories of chiral fermions belong to the most intrigu- 
ing properties of quantum field theory. Comprehensive reviews on anomalies can be 
found in the textbooks |jj 121 . 

Hydrodynamics is an ancient subject. Even in its relativistic form it appeared that 
everything relevant to its formulation could be found in Apart from stability is- 
sues that were addressed in the 1960s and 1970s 15] 15) leading to a second order 
formalism there seemed little room for new discoveries. The last years witnessed 
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however an unexpected and profound development of the formulation of relativistic 
hydrodynamics. The second order contributions have been put on a much more sys- 
tematic basis applying effective field theory reasoning |7 |. The lessons learned from 
applying the AdS/CFT correspondence [SJ to the plasma phase of strongly coupled 
non-abelian gauge theories ||9l \lU\ \TT\ played a major role (see lfT2l for a recent 
review). 

The presence of chiral anomalies in otherwise conserved currents has pro- 
found implications for the formulation of relativistic hydrodynamics. The trans- 
port processes related to anomalies have surfaced several times and independently 
lfT3l [141 [TSll . The axial current was the focus in |16| and the first application of 
the AdS/CFT correspondence to anomalous hydrodynamics can be found already 
in fT7 |. The full impact anomalies have on the formulation of relativistic hydrody- 
namics was however not fully appreciated until recently. 

The renewed interest in the formulation of relativistic hydrodynamics has its ori- 
gin mostly in the spectacular experimental evidence for collective flow phenomena 
taking place in the physics of heavy ion collisions at RHIC and LHC. These ex- 
periments indicate the creation of a deconfined quark gluon plasma in a strongly 
coupled regime. In the context of heavy ion collisions it was argued in [T8',T9 | that 
the excitation of topologically non-trivial gluon field configurations in the early non- 
equilibrium stages of a heavy ion collision might lead to an imbalance in the number 
of left- and right-handed quarks. This situation can be modeled by an axial chemical 
potential and it was shown that an external magnetic field leads to an electric current 
parallel to the magnetic field. This chiral magnetic effect leads then to a charge sep- 
aration perpendicular to the reaction plane in heavy ion collisions. The introduction 
of an axial chemical potential also allows to define a chiral magnetic conductivity 
which is simply the factor of proportionality between the magnetic field and the 
induced electric current. This effect is a direct consequence of the axial anomaly. 

The application of the fluid/gravity correspondence to theories including chiral 
anomalies lead to another surprise: it was found that not only a magnetic field in- 
duces a current but that also a vortex in the fluid leads to an induced current Il20ll2n . 
This is the chiral vortical effect. Again it is a consequence of the presence of a chi- 
ral anomaly. It was later realized that the chiral magnetic and vortical conductivities 
are almost completely fixed in the hydrodynamic framework by demanding the exis- 
tence of an entropy current with positive definite divergence |22 1. That this criterion 
did not fix the anomalous transport coefficients completely was noted in ll23l and 
various terms depending on the temperature instead of the chemical potentials were 
shown to be allowed as undetermined integration constants. See also |24| for a re- 
cent discussion of these anomaly coefficients with applications to heavy ion physics. 

In the meanwhile Kubo formulae for the chiral magnetic conductivity |25 1 and 
the chiral vortical conductivity |26| had been developed. Up to this point only 
pure gauge anomalies had been considered to be relevant since the mixed gauge- 
gravitational anomaly in four dimensions is of higher order in derivatives and was 
thought not to be able to contribute to hydrodynamics at first order in derivatives. 
Therefore it came as a surprise that in the application of the Kubo's formula for the 
chiral vortical conductivity to a system of free chiral fermions a purely temperature 
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dependent contribution was found. This contribution was consistent with some the 
earlier found integration constants and it was shown to arise if and only if the system 
of chiral fermions features a mixed gauge-gravitational anomaly 12711 . In fact these 
contributions had been found already very early on in lfT3l . The connection to the 
presence of anomalies was however not made explicit at that time. The gravitational 
anomaly contribution to the chiral vortical effect was also established in a strongly 
coupled AdS/CFT approach and precisely the same result as at weak coupling was 
found ESI. 

The argument based on a positive definite divergence of the entropy current 
allows to fix the contributions form pure gauge anomalies uniquely and provides 
therefore a non-renormalization theorem. No such result is known thus far for the 
contributions of the gauge-gravitational anomaly. 

A gas of weakly coupled Weyl fermions in arbitrary dimensions has been stud- 
ied in ||29l and confirmed that the anomalous conductivities can be obtained di- 
rectly from the anomaly polynomial under substitution of the field strength with 
the chemical potential and the first Pontryagin density by the negative of the tem- 
perature squared ll30l . Recently the anomalous conductivities have also been ob- 
tained in effective action approaches f3Tl l32l|. The contribution of the mixed gauge- 
gravitational anomaly appear on all these approaches as undetermined integrations 
constants. 

We will review here what can be learned from the calculation of the anomalous 
conductivities via Kubo formulae. The advantage of the usage of Kubo formulae is 
that they capture all contributions stemming either from pure gauge or from mixed 
gauge-gravitational anomalies. The disadvantage is that the calculations can be per- 
formed only with a particular model and only in a weak or in the gravity dual of the 
strong coupling regime. Along the way we will explain our point of view on some 
subtle issues concerning the definition of currents and of chemical potentials when 
anomalies are present. These subtleties lead indeed to some ambiguous results ll33l 
and |f34l . A first step to clarify these issues was done in |35| and a more general 
exposition of the relevant issues has appeared in |36|. 

The review is organized as follows. In section two we will briefly summarize the 
relevant issues concerning anomalies. We recall how vector like symmetries can al- 
ways be restored by adding suitable finite counterterms to the effective action ||39]| . 
A related but different issue is the fact that currents can be defined either as con- 
sistent or as covariant currents. The hydrodynamic constitutive relations depend on 
what definition of current is used. We discuss subtleties in the definition of the chem- 
ical potential in the presence of an anomaly and define our preferred scheme. We 
discuss the hydrodynamic constitutive relations and derive the Kubo formulae that 
allow the calculation of the anomalous transport coefficients from two point corre- 
lators of an underlying quantum field theory. 

In section three we apply the Kubo formulae to a theory of free Weyl fermions 
and show that two different contributions arise. They are clearly identifiable as being 
related to the presence of pure gauge and mixed gauge-gravitational anomalies. 

In section four we define a holographic model that implements the mixed gauge- 
gravitational anomaly via a mixed gauge-gravitational Chern-Simons term. We cal- 
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culate the same Kubo formulae as at weak coupling, obtaining the same values for 
chiral axi-magnetic and chiral vortical conductivities as in the weak coupling model. 

We conclude this review with some discussions and outlook to further develop- 
ments. 



2 Anomalies and Hydrodynamics 

In this section we will review briefly anomaUes. We compare the consistent with 
the covariant form of the anomaly and we introduce the Bardeen counterterm that 
allows to restore current conservation for vector Uke symmetries. Then we turn to 
the question of what we mean when we talk about the chemical potential. Two ways 
of introducing chemical potential, either through a deformation of the Hamiltonian 
or by demanding twisted boundary conditions along the thermal circle are shown 
to be in-equivalent in presence of an anomaly. Equivalence can still be achieved 
by introduction of a spurious axion field. We explain the implications for hologra- 
phy. The constitutive relations for anomalous currents are introduced in Landau and 
Laboratory frame. We discuss how they differ if we use the consistent instead of the 
covariant currents and derive the Kubo formulae for the anomalous conductivities. 



2.1 Anomalies 

Anomalies arise by integrating over chiral fermions in the path integral. They sig- 
nal a fundamental incompatibility between the symmetries present in the classical 
theory and the quantum theory. 

Unless otherwise stated we will always think of the symmetries as global sym- 
metries. But we still will introduce gauge fields. These gauge fields serve as classical 
sources coupled to the currents. As a side effect their presence promotes the global 
symmetry to a local gauge symmetry. It is still justified to think of it as a global 
symmetry as long as we do not introduce a kinetic Maxwell or Yang-Mills term in 
the action. 

In a theory with chiral fermions we define an effective action depending on these 
gauge fields by the path integral 



The vector field (x) couples to a classically conserved current —'Py^ . The 
charge operator Q can be the generator of a Lie group combined with chiral projec- 
tors ^± — 5 ( 1 ± 75 ) . General combinations are allowed although in the following 
we will mostly concentrate on a simple chiral U{1) symmetry for which we can take 
Q — 0^+. The fermions are minimally coupled to the gauge field and the classical 
action has an underUng gauge symmetry 




(1) 
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5f = -;A(x)ef , 5A^{x) ^ D^X{x) , (2) 

with denoting the gauge covariant derivative. Assuming that the theory has a 
classical limit the effective action in terms of the gauge fields allows for an expan- 
sion in h 

Wef / = Wo + nWi +Ti^W2 + ... (3) 

We find convenient to use the language of BRST symmetry by promoting the gauge 
parameter a ghost field c(x).[^The BRST symmetry is generated by 

sAn=DfiC , sc — —ic^ . (4) 

It is nilpotent = 0. The statement that the theory has an anomaly can now be 
neatly formalized. Since on gauge fields the BRST symmetry acts just as the gauge 
symmetry, gauge invariance translates into BRST invariance. An anomaly is present 
if 

sWeff = and £/ ^ sB . (5) 

Because of the nilpotency of the BRST operator the anomaly has to fulfill the Wess- 
Zumino consistency condition 

s£/^0. (6) 

As indicated in (j5]l this has a possible trivial solution if there exists a local functional 
B[A^] such that sB — si/. An anomaly is present if no such B exists. The anomaly 
is a quantum effect. If it is of order h" and if a suitable local functional B exists we 
could simply redefine the effective action as W^ff = W^ff — B and the new effec- 
tive action would be BRST and therefore gauge invariant. The form and even the 
necessity to introduce a functional B might depend on the particular regularization 
scheme chosen. As we will explain in the case of an axial and vector symmetry a 
suitable B can be found that always allows to restore the vectorlike symmetry, this 
is the so-called Bardeen counterterm |39|. The necessity to introduce the Bardeen 
counterterm relies however on the regularization scheme chosen. In schemes that au- 
tomatically preserve vectorlike symmetries, such as dimensional regularization, the 
vector symmetries are automatically preserved and no counterterm has to be added. 
Furthermore the Adler-Bardeen theorem guarantees that chiral anomalies appear 
only at order fi. Their presence can therefore by detected in one loop diagrams such 
as the triangle diagram of three currents. 

We have introduced the gauge fields as sources for the currents 

Weff[A] = {J^). (7) 



5A^{x) 



For chiral fermions transforming under a general Lie group generated by T" the 
chiral anomaly takes the form |JJ 



A recent comprehensive review on BRST symmetry is 1401 
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T"d^ I AydpAfj + —AyApAfj 



(8) 



Where 7] = +1 for left-handed fermions and Tj = — 1 for right-handed fermions. 
Differentiating with respect to the ghost field (the gauge parameter) we can derive a 
local form. To simplify the formulas we specialize this to the case of a single chiral 
U{1) symmetry taking T" = 1 



pa ■ 



(9) 



This is to be understood as an operator equation. Sandwitching it between the vac- 
uum state |0) and further differentiating with respect to the gauge fields we can 
generate the famous triangle form of the anomaly 



(10) 



The one point function of the divergence of the current is non-conserved only in the 
background of parallel electric and magnetic fields whereas the non-conservation of 
the current as an operator becomes apparent in the triangle diagram even in vacuum. 

By construction this form of the anomaly fulfills the Wess-Zumino consistency 
condition and is therefore called the consistent anomaly. In analogy we call the 
current defined by (|7| the consistent current. 

For a f/ ( 1 ) symmetry the functional differentiation with respect to the gauge field 
and the BRST operator s commute. 



5Af, (x) 



= 0. 



(11) 



An immediate consequence is that the consistent current is not BRST invariant but 
rather obeys 



sJ^ 



(12) 



where we introduced the Chern-Simons current — e'^^P'^AyPpi with d^K^ — 
yf^'P'-F^.Fp^. 

With the help of the Chern-Simons current it is possible to define the so-called 
covariant current (in the case of a symmetry rather the invariant current) 



1 



24;r2 



fulfilhng 



sP^Q. 

The divergence of the covariant current defines the covariant anomaly 



(13) 
(14) 
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V = ^e^^''%v^pa. (15) 

Notice that the Chem-Simons current cannot be obtained as the variation with re- 
spect to the gauge field of any funtional. It is therefore not possible to define an 
effective action whose derivation with respect to the gauge field gives the covariant 
current. 

Let us suppose now that we have one left-handed and one right-handed fermion 
with the corresponding left- and right-handed anomalies. Instead of the left-right ba- 
sis it is more convenient to introduce a vector-axial basis by defining the vectorlike 
current Jy — +Jg and the axial current = ^ ^r - Let ^ni^) be the gauge 
field that couples to the vectorlike current and A|j (x) be the gauge field coupling to 
the axial current. The (consistent) anomalies for the vector and axial current turn 
out to be 

« = ^2^''''iF>pX+F',vF^x)- (17) 

As long as the vectorlike current corresponds to a global symmetry nothing has 
gone wrong so far If we want to identify the vectorlike current with the electric- 
magnetic current in nature we need to couple it to a dynamical photon gauge field 
and now the non-conservation of the vector current is worrysome to say the least. 
The problem arises because implicitly we presumed a regularization scheme that 
treats left-handed and right-handed fermions on the same footing. As pointed out 
first by Bardeen this flaw can be repaired by adding a finite counterterm (of order h) 
to the effective action. This is the so-called Bardeen counterterm and has the form 

Adding this counterterm to the effective action gives additional contributions of 
Chern-Simons form to the consistent vector and axial currents. With the particular 
coefficient chosen it turns out that the anomaly in the vector current is canceled 
whereas the axial current picks up an additional contribution such that after adding 
the Bardeen counterterm the anomalies are 

Vy=0, (19) 

d.J'A = ^2<^''"H^f;:,F^x +F^vF^x) • (20) 

This definition of currents is mandatory if we want to identify the vector current with 
the usual electromagnetic current in nature! It is furthermore worth to point out that 
both currents are now invariant under the vectorlike U{1) symmetry. The currents 
are not invariant under axial transformation, but these are anomalous anyway. 

Generalizations of the covariant anomaly and the Bardeen counterterm to the 
non-abelian case can be found e.g. in |[T|. 
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There is one more anomaly that will play a major role in this review, the mixed 
gauge-gravitational anomaly ["^Tl.l^So far we have considered only spin one cur- 
rents and have coupled them to gauge fields. Now we also want to introduce the 
energy-momentum tensor through its coupling to a fiducial background metric g^v 
Just as the gauge fields, the metric serves primarily as the source for insertions of 
the energy momentum tensor in correlation functions. Just as in the case of vector 
and axial currents, the mixed gauge-gravitational anomaly is the statement that it is 
impossible in the quantum theory to preserve at the same time the vanishing of the 
divergence of the energy-momentum tensor and of chiral (or axial) U{1) currents. It 
is however possible to add Bardeen counterterms to shift the anomaly always in the 
sector of the spin one currents and preserve translational (or diffeomorphism) sym- 
metry. If we have a set of left-handed and right-handed chiral fermions transform- 
ing under a Lie Group generated by (r„)i and (r^)^ in the background of arbitrary 
gauge fields and metric, the anomaly equations in covariant form are 

Df^T^^^O, (21) 

The purely group theoretic factors are 

dabc = ^tr(r„{7;, r,})i - ^ti{Ta{Tb, T,})r , (23) 
ba^tr{Ta)L-tr{Ta)R. (24) 

Anomalies are completely absent if and only if dabc — and ba=Q. 



2.2 Chemical Potentials for anomalous symmetries 

Thermodynamics of systems with conserved charges can be described in a grand 
canonical ensemble where a Lagrange multiplier jj, ensures that the partition func- 
tion fulfills 

r^He>. (25) 

The textbook approach is to consider a deformation of the Hamiltonian 

H^H-iiQ, (26) 

where Q is the charge in question. We can think of this as arising from the coupling 
of the (fiducial) gauge field to the current y and giving a vacuum expectation 
value to Aq = jj.. Since the fiducial gauge field leads to local gauge invariance we 



^ In D = 4k + 2 dimensions also purely gravitational anomalies can appear 1421 . 
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can remove the coupling in the Hamiltonian by the gauge transformation Aq 
Aq + doX with X = -pit. 



Fig. 1 At finite temperature field theories are defined on the Keldysh-Schwinger contour in the 
complexified time plane. The initial stated at f, is specified through the boundary conditions on the 
fields. The endpoint of the contour is at /, — ;/3 where ji = l/T. 



In the context of finite temperature field theory such a gauge transformation is 
however not really allowed. One needs to define the field theory on the Keldysh- 
Schwinger contour in the complexified time plane as shown in figure ([TJ. Fields are 
taken to be periodic or anti-periodic along the imaginary time direction t = — /t with 
period P — l/T where T is the temperature 

f(r,-/i3) = ±f(f,), (27) 

with the plus sign for bosons and the minus sign for fermions. The gauge transfor- 
mation that removes the constant zero component of the gauge field is not periodic 
along the contour and therefore changes the boundary conditions on the fields. After 
the gauge transformation with x = —pi-t the fields obey the boundary conditions 

f (f,-/j3) =±£''"''^'f^(f,). (28) 

Demanding these "twisted" boundary conditions is of course completely equiva- 
lent to having Aq = /X . The gauge invariant statement is that a charged field parallel 
transported around the Keldysh-Schwinger contour picks up a factor of exp(g'/ij3). 
As long as we have honest non-anomalous symmetries under consideration we have 
therefore two (gauge)-equivalent formaUsms of how to introduce the chemical po- 
tential summarized in table [T] ll37ll . 



Table 1 Two formalisms for the chemical potential 



Formalism 


Hamiltonian 


Boundary condition 


(A) 




f(f,-'|8) = ±'f'(f,) 


(B) 


H 
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One convenient point of view on formalism (B) is the following. In a real time 
Keldysh-Schwinger setup we demand some initial conditions at initial (real) time 
t — ti. These initial conditions are given by the boundary conditions in (B). From 
then on we do the (real) time development with the microscopic Hamiltonian H. In 
principle there is no need for the Hamiltonian H to preserve the symmetry present at 
times t <ti. This seems an especially suited approach to situations where the charge 
in question is not conserved by the real time dynamics. In the case of an anomalous 
symmetry we can start at f = f,- with a state of certain charge. As long as we have 
only external gauge fields present the one-point function of the divergence of the 
current vanishes and the charge is conserved. This is not true on the full theory 
since even in vacuum the three-point correlators are sensitive to the anomaly. For 
the formulation of hydrodynamics in external fields the condition that the one-point 
functions of the currents are conserved as long as there are no parallel electric and 
magnetic external fields (or a metric that has non-vanishing Pontryagin density) is 
sufficient. 

Let us assume now that Q is an anomalous charge, i.e. its associated current 
suffers from chiral anomalies. We first consider formalism (B) and ask what happens 
if we do now the gauge transformation that would bring us to formalism (A). Since 
the symmetry is anomalous the action transforms as 

S[A + dx]=S[A]+ J d^xxe^'^'P^ {CiF^yFp^+C2R"p^vR^apx) , (29) 

with the anomaly coefficients Ci and Ci depending on the chiral fermion content. It 
follows that formalisms (A) and (B) are physically inequivalent now, because of the 
anomaly. However, we would like to still come as close as possible to the formalism 
of (A) but in a form that is physically equivalent to the formalism (B). To achieve 
this we proceed by introducing a non-dynamical axion field (x) and the vertex 

Se[A,0] = J d^xQe^^'P^ {c.F^vFpx +C2R" p^yR^ apx) ■ (30) 

If we demand now that the "axion" transforms as© ^ © — x under gauge transfor- 
mations we see that the action 

S,ot[A,©]=S[A]+S@[A,©] (31) 

is gauge invariant. Note that this does not mean that the theory is not anomalous 
now. We introduce it solely for the purpose to make clear how the action has to 
be modified such that two field configurations related by a gauge transformation 
are physically equivalent. It is better to consider as coupling and not a field, i.e. 
we consider it a spurion field. The gauge field configuration that corresponds to 

' If dynamical gauge fields are present, such as the gluon fields in QCD even the one point function 
of the charge does decay over (real) time due to non-perturbative processes (instantons) or at finite 
temperature due to thermal sphaleron processes 1 38 1. Even in this case in the limit of large number 
of colors these processes are suppressed and can e.g. not be seen in holographic models in the 
supergravity approximation. 



Anomalous Transport from Kubo Formulae 1 1 

formalism (B) is simply Aq = 0. A gauge transformation with x = l^t on the gauge 
invariant action Star makes clear that a physically equivalent theory is obtained by 
chosing the field configuration Aq = ji and the time dependent coupling ~ —jxt. 
If we define the current through the variation of the action with respect to the gauge 
field we get an additional contribution from Sq, 

4=4Cie^''P^dy0Fp^, (32) 

and evaluating this for & — —jit we get the spatial current 

J'g=4CiliB'". (33) 

We do not consider this to be the chiral magnetic effect! This is only the contribution 
to the current that comes from the new coupling that we are forced to introduce by 
going to formalism (A) from (B) in a (gauge)-equivalent way. As we will see in 
the following chapters the chiral magnetic and vortical effect are on the contrary 
non-trivial results of dynamical one-loop calculations. 

What is the Hamiltonian now based on the modified formalism (A)? We have to 
take of course the new coupling generated by the non-zero 0. The Hamiltonian now 
is therefore 

H-iJ.(^Q + 4Ci J d\e^'j%djAk^ , (34) 

where for simplicity we have ignored the contributions from the metric terms. 

For explicit computations in sections [3] and |4] we will introduce the chemical 
potential through the formalism (B) by demanding twisted boundary conditions. It 
seems the most natural choice since the dynamics is described by the microscopic 
Hamiltonian H. The modified (A) based on the Hamiltonian ( |34j i is however not 
without merits. It is convenient in holography where it allows vanishing temporal 
gauge field on the black hole horizon and therefore a non-singular Euclidean black 
hole geometry.]^ 



2.2.1 Hydrodynamics and Kubo formulae 

The modern understanding of hydrodynamics is as an effective field theory. The 
equations of motion are the (anomalous) conservation laws of the energy-momentum 
tensor and spin one currents. These are supplemented by expression for the energy- 
momentum tensor and the current which are organized in a derivative expansion, 
the so-called constitutive relations. Symmetries constrain the possible terms. In the 
presence of chiral anomalies the constitutive relations for the energy-momentum 
tensor and the currents in the Landau frame are 



It is possible to define a generalized formalism to make any choice for the gauge field Aq = v, 
so that one recovers formalism (A) when v = ^ and formalism (B) when v = as particular cases 
(see 1431 for details). 
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= eu^u^ + pP^" - riP^^P^^aaii - ^P^^dau" , (35) 
J^^p,u^+L,,„[Ej;-TP^'^Do,^)+CB^,+^:cor (36) 

It is important to specify that these are the constitutive relations for the covariant 
currents! Here e is the energy density, p the pressure density, the local fluid veloc- 
ity. P^"^ = g^^ +u^u^ is the transverse projector to the fluid velocity. (7^v is the sym- 
metric traceless shear tensor The non-anomalous transport coefficients are the shear 
viscosity T], the bulk viscosity ^ and the electric conductivities External electric 
and magnetic fields are covariantized via = F^^ and Ba = ^e^^P^UyF„ p^. 

The anomalous transport coefficients are the chiral magnetic conductivities (^J^ 
and the chiral vortical conductivities . At first order in derivatives the notion of 
fluid velocity is ambiguous and needs to be fixed by prescribing a choice of frame. 
We remark that the constitutive relations ( (35| ) and ( [36) 1 are valid in the Landau frame 
where T^^uy ~ eu^. 

To compute the Kubo formulae for the anomalous transport coefficients it turns 
out that the Landau frame is not the most convenient one. It fixes the definition 
of the fluid velocity through energy transport. Transport phenomena related to the 
generation of an energy current are therefore not directly visible, rather they are 
absorbed in the definition of the fluid velocity. It is therefore more convenient to go 
to another frame in which we demand that the definition of the fluid velocity is not 
influenced when switching on an external magnetic field or having a vortex in the 
fluid. In such a frame the constitutive relations take the form 

T^^ = {e + p)u^u'' ^ pg^" - riP^^P"^ Oafi - l^P^" dau"" + Q^u'' + Q''u^ ,{31) 
= a§B^ + a^co^ , (38) 

^ pu^+E(E^-TPf'"Da(^j:))+aBB^ + av(0^. (39) 

In order to avoid unnecessary clutter in the equations we have specialized now to a 
single U{1) charge. Notice that now there is a sort of "heat" current present in the 
constitutive relation for the energy-momentum tensor 

The derivation of Kubo formulae is better based on the usage of the consistent 
currents. Since the covariant and consistent currents are related by adding suitable 
Chern-Simons currents the constitutive relations for the consistent current receives 
additional contribution from the Chern-Simons current 

J^^J^--^Kr (40) 

If we were to introduce the chemical potential according to formalism (A) via a 
background for the temporal gauge field we would get an additional contribution to 
the consistent current from the Chern-Simons current. In this case it is better to go to 
the modified formalism (A') that also introduces a spurious axion field and another 
contribution to the current J@ ((33]l has to be added 
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J'=J'-^2K'+4- (41) 

For the derivation of the Kubo formulae it is therefore more convenient to work 
with formalism (B) in which Aq = and the chemical potential is introduced via the 



boundary conditions (28 1. Otherwise there arise additional contributions to the two 
point functions. We will briefly discuss them in the next subsection. 

From the microscopic view the constitutive relations should be interpreted as the 
one-point functions of the operators and in a near equilibrium situation, i.e. 
gradients in the fluid velocity, the temperature or the chemical potentials are as- 
sumed to be small. From this point of view Kubo formulae can be derived. In the 
microscopic theory the one-point function of an operator near equilibrium is given 
by linear response theory whose basic ingredient are the retarded two-point func- 
tions. If we consider a situation with only an external electric field in z-direction 
and all other sources switched off, i.e. the fluid being at rest = (1,0,0,0) and no 
gradients in temperature or chemical potentials the constitutive relations are simpli- 
fied to 

r = EE' . (42) 

The electric field is = icoA- in terms of the vector potential and using linear 
response theory the induced current is given through the retarded two-point function 
by 

f'^{J-J')A''. (43) 

Equating the two expressions for the current we find the Kubo's formula for the 
electric conductivity 

i: = lim — If'P) . (44) 

ffl^O CO 

This has to be evaluated at zero momentum. The limit in the frequency follows 
because the constitutive relation are supposed to be valid only to lowest order in the 
derivative expansion, therefore one needs to isolate the first non-trivial term. 

Now we want to find some simple special cases that allow the derivation of 
Kubo's formulae for the anomalous conductivities. A very convenient choice is to 
go to the restframe = (1,0,0,0), switch on a vector potential in the y-direction 
that depends only on the z direction and at the same time a metric deformation 
gjiv = Tl/iv + h/xy with the only non-vanishing component hoy depending on z 
only. To linear order in the background fields the non-vanishing components of the 
energy-momentum tensor and the current are 

= -aldAy - a^dJiQy , (45) 

r = -GBd^Ay - Gvd-JlQy , (46) 

since in the formalism {B) neither the Chern-Simons term nor the coupling con- 
tribute. Going to momentum space and differentiating with respect to the sources Ay 
and /iQy we find therefore the Kubo formulae lfT9ll26l 
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Ob = limA.._ 








C7| = limk__ 









(47) 



All these correlators are to be taken at precisely zero frequency. As these for- 
mulas are based on linear response theory the correlators should be understood as 
retarded ones. They have to be evaluated however at zero frequency and therefore 
the order of the operators can be reversed. From this it follows that the chiral vorti- 
cal conductivity coincides with the chiral magnetic conductivity for the energy flux 

We also want to discuss how these transport coefficients are related to the ones 
in the more commonly used Landau frame. They are connected by a redefinition of 
the fluid velocity of the form 



e+p 



(48) 



to go from (37i-(39i to (35 i-([36]l. The corresponding transport coefficients of the 
Landau frame are therefore 



Ib = lim — ^re„H ( (// , 

k„^o2k„f;f \\ / £ + p 



k.l 



e + p 



(49) 
(50) 



where we have employed a slightly more covariant notation. The generalization to 
the non-abelian case is straightforward. 

It is also worth to compare to the Kubo formulae for the dissipative transport 



coefficients as the electric conductivity ( 44 1. In the dissipative cases one first goes 
to zero momentum and then takes the zero frequency limit. In the anomalous con- 
ductivities this is the other way around, one first goes to zero frequency and then 
takes the zero momentum limit. Another observation is that the dissipative transport 
coefficients sit in the anti-Hermitean part of the retarded correlators, i.e. the spectral 
function whereas the anomalous conductivities sit in the Hermitean part. The rate at 
which an external source // does work on a system is given in terms of the spectral 
function of the operator coupling to that source as 



dW 



1 



cof,{-co)p"{co)fj{co). 



(51) 



Notice that hoy can also be understood as the so-called gravito-magnetic vctor potential A^, which 
is related to the gravito-magntic field by = V x Ag. This allows to interpret Oy not only as the 
generation of a current due to a vortex in the fluid, i.e. the chiral vortical effect, but also as a chiral 
gravito-magnetic conductivity giving rise to a chiral gravito-magnetic effect, see |44| for details. 



Anomalous Transport from Kubo Formulae 



15 



The anomalous transport phenomena therefore do no work on the system, first they 
take place at zero frequency and second they are not contained in the spectral func- 
tion p = ^ (G, —Gl). 



2.3 Contributions to the Kubo formulae 

Now we want to give a detailed analysis of the different Feynman graphs that con- 
tribute to the Kubo formulae in the different formalisms for the chemical potentials. 
The simplest and most economic formalism is certainly the one labeled (B) in which 
we introduce the chemical potentials via twisted boundary conditions. The Hamilto- 
nian is simply the microscopic Hamiltonian H. Relevant contributions arise only at 
first order in the momentum and at zero frequency and in this kinematic limit only 
the Kubo formulae for the chiral magnetic conductivity is affected. In the figure Q 
we summarize the different contributions to the Kubo's formulas in the three ways 
to introduce the chemical potential. 




finite T. // loop vacuum loop coupling to spurious 9 Beld 

Fig. 2 Contributions to the Kubo's formula for the chiral magnetic conductivity in the different 
formahsms for the chemical potential. 

The first of the Feynman graphs is the same in all formalisms. It is the genuine 
finite temperature and finite density one-loop contribution. This graph is finite be- 
cause the Fermi-Dirac distributions cutoff the UV momentum modes in the loop. In 
the formalism (A) we need to take into account that there is also a contribution from 
the triangle graph with the fermions going around the loop in vacuum, i.e. without 
the Fermi-Dirac distributions in the loop integrals. For a non-anomalous symmetry 
this graph vanishes simply because on the upper vertex of the triangle sits a field 
configuration that is a pure gauge. If the symmetry under consideration is however 
anomalous the triangle diagram picks up just the anomaly. Even pure gauge field 
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configurations become physically distinct from the vacuum and therefore this dia- 
gram gives a non-trivial contribution. On the level of the constitutive relations this 



contribution corresponds to the Chern-Simons current in (40 1. We consider this con- 
tribution to be unwanted. After all the anomaly would make even a constant value 
of the temporal gauge field Aq observable in vacuum. An example is provided for a 
putative axial gauge field A^. If present the absolute value of its temporal compo- 
nent would be observable through the axial anomaly. We can be sure that in nature 
no such background field is present. The third line (A') introduces also the spuri- 
ous axion field the only purpose of this field is to cancel the contribution from 
the triangle graph. This cancellation takes place by construction since (A') is gauge 
equivalent to (B) in which only the first genuine finite T,ii part contributes. It cor- 



responds to the contribution of the current Jq in (141 



We further emphasize that these considerations are based on the usage of the 
consistent currents. In the interplay between axial and vector currents additional 
contributions arise from the Bardeen counterterm. It turns out that the triangle or 
Chern-Simons current contribution to the consistent vector current in the formal- 
ism (A) cancels precisely the first one lf34l [35l . Our take on this is that a constant 
temporal component of the axial gauge field Aq = /ij would be observable in nature 
and can therefore be assumed to be absent. The correct way of evaluating the Kubo 
formulae for the chiral magnetic effect is therefore the formalism (B) or the gauge 
equivalent one (A'). 

At this point the reader might wonder why we introduced yet another formal- 
ism (A') which achieves appearently nothing but being equivalent to formalism (B). 
At least from the perspective of holography there is a good reason for doing so. 
In Holography the strong coupling duals of gauge theories at finite temperature in 
the plasma phase are represented by five dimensional asymptotically Anti- de Sit- 
ter black holes. Finite charge density translates to charged black holes. These black 
holes have some non-trivial gauge flux along the holographic direction represented 
by a temporal gauge field configuration of the form Ao(r) where r is the fifth, holo- 
graphic dimension. It is often claimed that for consistency reasons the gauge field 
has to vanish on the horizon of the black hole and that its value on the boundary can 
be identified with the chemical potential 

Ao(r//)=0 and Ao(r^oc)=^. (52) 

According to the usual holographic dictionary the gauge field values on the bound- 
ary correspond to the sources for currents. A non-vanishing value of the temporal 
component of the gauge field at the boundary is therefore dual to a coupling that 



modifies the Hamiltonian of the theory just as in (26 1. Thus with the boundary con 



ditions (52i we have the holographic dual of the formalism (A). If anomalies are 
present they are represented in the holographic dual by five-dimensional Chern- 
Simons terms of the form A /\F AF. The two point correlator of the (consistent) 
currents receives now contributions from the Chern-Simons term that is precisely of 
the form of the second graph in (A) in figure |2] As we have argued this is an a priory 
unwanted contribution. We can however cure that by introducing an additional term 
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in the action of the form ( [30| l Hving only on the boundary of the holographic space- 
time. In this way we can implement the formalism (A'), cancel the unwanted triangle 
contribution with the third graph in (A') in figure|2]and maintain Ao(r//) = 0! 

The claim that the temporal component of the gauge field has to vanish at the 
horizon is of course not unsubstantiated. The reasoning goes as follows. The Eu- 
clidean section of the black-hole space time has the topology of a disc in the r, T 
directions, where T is the Euclidean time. This is a periodic variable with period 
P = 1/T where T is the (Hawking) temperature of the black hole and at the same 
time the temperature in the dual field theory. Using Stoke's law we have 



AodT^ / FrodrdT, (53) 

do Jd 

where /vo is the electric field strength in the holographic direction and D is a Disc 
with origin at r = r// reaching out to some finite value of rf. If we shrink this disc 
to zero size, i.e. let rf — > the r.h.s. of the last equation vanishes and so must the 
I.h.s. which approaches the value j3Ao(r//). This implies that Ao(r//) = 0. If on the 
other hand we assume that Ao(r//) 7^ then the field strength must have a delta type 
singularity there in order to satisfy Stokes theorem. Strictly speaking the topology 
of the Euclidean section of the black hole is not anymore that of a disc since now 
there is a puncture at the horizon. It is therefore more appropriate to think of this as 
having the topology of a cylinder Now if we want to implement the formalism (B) 
in holography we would find the boundary conditions 

Ao(r//)=^ and Ao(r^oo)=0, (54) 

and precisely such a singularity at the horizon would arise. In addition we would 
need to impose twisted boundary conditions around the Euclidean time T for the 



fields just as in (28 1. Now the presence of the singularity seems to be a good thing: 
if the space time would still be smooth at the horizon it would be impossible to 
demand these twisted boundary conditions since the circle in T shrinks to zero size 
there. If this is however a singular point of the geometry we can not really shrink 
the circle to zero size. The topology being rather a cylinder than a disc allows now 
for the presence of the twisted boundary conditions. 

It is also important to note that in all formalisms the potential difference between 
the boundary and the horizon is given by fi. This has a very nice intuitive interpre- 
tation. If we bring a unit test charge from the boundary to the horizon we need the 
energy AE = jx. In the dual field theory this is just the energy cost of adding one unit 
of charge to the thermalized system and coincides with the elementary definition of 
the chemical potential. 

From now on we will always only consider the genuine finite T, /i contribution 
that is the only one that arises in formalism (B). 

The rest of this review is devoted to the explicit evaluation of these Kubo's for- 
mulae in two different systems: free chiral fermions and a holographic model imple- 
menting the chiral and gravitational anomalies by suitable five dimensional Chem- 
Simons terms. 



18 Karl Landsteiner, Eugenio Megias and Francisco Pena-Benitez 




Horizon Boundary 



Fig. 3 A sketch of the Euclidean black hole topology. A singularity at the horizon arises if we do 
not choose the temporal component of the gauge field to vanish there. On the other hand allowing 
the singularity to be present changes the topology to the one of a cylinder and this in turn allows 
twisted boundary conditions. 



3 Weyl fermions 

We will now evaluate the Kubo formulae for the chiral magnetic, chiral vortical and 



energy flux conductivities (47 i for a theory of free chiral fermions f transform- 
ing under a global symmetry group G generated by matrices {TaY g. 

We denote the generators in the Cartan subalgebra by Chemical potentials pia 
can be switched on only in the Cartan subalgebra. Furthermore the presence of the 
chemical potentials breaks the group G to a subgroup G. Only the currents that lie 
in the unbroken subgroup are conserved (up to anomalies) and participate in the hy- 
drodynamics. The chemical potential for the fermion is given by /i^ = lal^a, 
where we write the Cartan generator //„ ~ qliS^ g in terms of its eigenvalues, the 
charges q 
fulfilUng 



f ^ f 

charges qi, . The unbroken symmetry group G is generated by those matrices Tj g 



Tigii^ = llfTig. (55) 

There is no summation over indices in the last expression. From now on we will 
assume that all currents !„ lie in directions indicated in ( |55] l. We define the chemical 
potential through the boundary condition on the fermion fields around the thermal 
circle, i.e. we adopt the formalism (B) discussed in previous section, 

f^(T-i8) = -e'5''V-^(T). (56) 

Therefore the eigenvalues of are /(»„ + /i^ for the fermion spiecies / with 
d),, = nT{2n+ 1) the fermionic Matsubara frequencies. A convenient way of ex- 
pressing the current and the energy-momentum tensor is in terms of Dirac fermions 
and writing 

J'a^ t T!f%Y^+^f , = ^ f %{fd' + Yd^)3^+Wf , (57) 

/..?=1 /=1 



where we used the chiral projector ^± = ^ (1 ± 75). The fermion propagator is 
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q + k 




%/g = ^L4,(/a)^,q)^+7^^f , Micbf.q)^ (58) 

with io)^ — i&)„ + ji^, = {l,tq), q — J- and Eg — \q\. For simplicity in the ex- 
pressions we consider only left-handed fermions, but one can easily include right- 
handed fermions as well as they contribute in all our calculations in the same way 
as the left-handed ones up to a relative minus sign. 

We will address in detail the computation of the vortical conductivities and sketch 
only the calculation of the magnetic conductivities since the latter one is a trivial 
extension of the calculation of the chiral mangetic conductivity in ll25l . Then we 
show the results for the other conductivity coefficients. 



3.1 Chiral Vortical Conductivity 



The vortical conductivity is defined from the retarded correlation function of the 
current and the energy momentum tensor or energy current r'''(x') (57 1, i.e. 



1 

2 ^'J" ' 



(59) 



Going to Fourier space, one can evaluate this quantity as 



/=1 



wf 



d q 

W) 



3 '^ijn 



Sf f(g)Y Sf f(q+k) [-fq'+fiStf) 



(60) 



which corresponds to the one loop diagram of figure]?] The vertex of the two quarks 
with the graviton is ^ 5^ g, and therefore we find only contributions from the diag- 
onal part of the group G. The metric we use through this section is the usual one 
in field theory computations, g^y — diag(l, — 1, — 1, — 1). We can split Gjj^ into two 
contributions, i.e. 

Grw-<(o,)(^)+Gr,ao)W> (61) 
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which correspond to the terms "f^q^ and yHS)^ in eq. ( 60 1 respectively. We will focus 
first on the computation of G^^q^.j. After computation of the Dirac trace in eq. (60 1, 
this term writes 



tM=± 



-q 



-'q+k 



tu 



A, (iW-^, k)ziu {io}^ + ifOn , q + k) . (62) 



At this point one can make a few simplifications. Note that due to the antisymmetric 



tensor e,y„, the two terms proportional to q' inside the bracket in eq. (62i vanish 



Regarding the term eij„q-^k', it leads to a contribution ^ eij„k-'k' after integration in 
d^q, which is zero. Then the only term which remains is the one not involving eij„. 
We can now perform the sum over fermionic Matsubara frequencies. One has 



^^40«-^,q)4„(/d)-^ + /a)„,q + p) = 

s>f 



tn(E^-tllf)-im{E^^l,-uflf)+^iu-t) 



(63) 



where n{x) — l/{e^^ + 1) is the Fermi-Dirac distribution function. In eq. (63 i we 
have considered that (On = 2nTn is a bosonic Matsubara frequency. This result is 
also obtained in Ref ||251 . After doing the analytic continuation, which amounts to 



replacing /(»„ by + '£ in eq. (63 1, one gets 



(k) 



a f 



E 

r.u=± 



d^q q^k„~{q-k)q„ 
un{Ec, - 



EqEqj^li 



' t^f) - tn{Ey+k - uji f) + 5 (f - m) 



ko + ie + tEa — uE, 



q+k 



(64) 



The term proportional to ^ j (f — m) corresponds to the vacuum contribution, and it 
is ultraviolet divergent. By removing this term the finite temperature and chemical 
potential behavior is not affected, and the result becomes ultraviolet finite because 
the Fermi-Dirac distribution function exponentially suppresses high momenta. By 
making both the change of variable q — ^ — q — k and the interchange m — > — f and 
t — M in the part of the integrand involving the term —tn{Eg^j^ — uii^), one can 



express the vacuum substracted contribution of eq. (64 1 as 



E " 

t,u=± 



A 1 ( 2 (q-k) 

(27r)3 EqEq+k V k2 
n{E„-llf)+n{Eq+pLf) 



ko + i£ + tEa + 



q+k 



(65) 



Anomalous Transport from Kubo Formulae 



21 



where we have used that n{Eg — tji^ ) +n{Eg+tjif) — n{Eq — jx-f )+n{Eq + /i-^) since 
t = ±1. The result has to be proportional to k„, so to reach this expression we have 
replaced q„ by (q • k)A;„/k^ in eq. (64 1. At this point one can perform the sum over 
Mbyusing^„^-i-M/(ai+Ma2) 



2fl2/ (a J - 



- 02), and the integration over angles by 



considering q • k = EqEj^x and E^^^. ^ E^ +E^ + lEqE^x, where x : = cos ( ) and 
is the angle between q and k. Then one gets the final result 



1 + ^ L N " + + tko)] log 

where 12, = ^0 + + tEq . and 

f{q)=n{Eq-pLf)+n{Eq + pLf). 

The steps to compute 



(66) 



^2^ 

£1} 



-kf 



(67) 



in eq. (61 



are similar. In this case the Dirac trace leads 
1 the only contribution comes from the trace 



to a different tensor structure, in whic : 
involving 75. The sum over fermionic Matsubara frequencies involves an extra iS)^ , 
i.e. 



^ ^ im^A, [imf ,q)Ai, [im^ + i(o„ , q + k) 



1 



I (Or, 



- tEq - uEqj^k 



(68) 



Eqn{Eq - t^l f) - {Eq+k " ui(0„)n{E q+j, - M/i' 



[Eq - Eq+k + Ui(0n 



The last term inside the bracket in the r.h.s. of eq. ( [68] ) corresponds to the vacuum 
contribution which we choose to remove, as it leads to an ultraviolet divergent con- 
tribution after integration in cfiq. Making similar steps as for G^^^q^j, i.e. performing 
the sum over u and integrating over angles, one gets the final result 



G 



am 



(k) 



i k. 



32n^ 



t=± 



Atqkh) - [k^ - kl) {2q + tko) log 



Q^-iq-ky 



where 



f^{qM)^qf{q)+tk)n{Eq + t^f). 



, (69) 



(70) 



The result for the vacuum substracted contribution of the retarded correlation func- 
tion of the current and the energy momentum tensor, G^{k), writes as a sum of 



eqs. ( [66] l and (69 1, according to eq. (61 1. From these expressions one can compute 
the zero frequency, zero momentum, limit. Since 
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lim lim > log — — :r- 



2k 

q ' 



the relevant integrals are 

dqqf{q) 



i: 



dqf^{q,kQ=0) 



in')' 



(71) 



(72) 



Finally it follows from eqs. (66 1 and (69 1 that the zero frequency, zero momentum, 
vortical conductivity writes 



{(yv)a 



1 

87r2 



E Tif 



2n^ 



£ti- (r„{//,„//,})Ali,M. + ^r^tr [To) 



(73) 



Both G^jQ^jand G^(jo) ^^^'^ '^^ '^^^ same contribution in [Ov) 



Eq. ( jTij was first 

derived in Il27ll . and it constitutes our main result in this section, . The term involving 
the chemical potentials is induced by the chiral anomaly. More interesting is the 
term ^ which is proportional to the gravitational anomaly coefficient ba [4Tll42]| . 
This means that a non-zero value of this term has to be attributed to the presence 
of a gravitational anomaly. The Matsubara frequencies &>„ = nT{2n+ 1) generate 
a dependence on nT in the final result as compared to the chemical potentials, and 



then no factors of n show up for the term ^ T in eq. (73 1. Right-handed fermions 



contribute in the same way but with a relative minus sign. Therefore the ~ T term 



appears only when the current in eq. ( 59 1 has an axial component. The correlator 
with a vector current does not have this gravitational anomaly contribution. 



3.2 Chiral Magnetic Conductivity 



The chiral magnetic conductivity in the case of a vector and an axial U{1) symmetry 
was computed at weak coupling in ll25l . The corresponding Kubo formula involves 



the two point function of the current, see first expression in eq. (47 1. Following 
the same method, we have computed it for the unbroken (non-abelian) symmetry 
group G. The relevant Green function is (TT\ 



P mf 



d q 

W) 



Sf f{q)YSf f{q + k)Y^ 



(74) 



The evaluation of this expression is exactly as in ||251 so we skip the details. The 
zero frequency, zero momentum, limit of the magnetic conductivity is 
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f.s=i 



Y,ti{T,{Tb,H,})^c. 
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(75) 



In the second equality of eq. (75 1 we have made use of eq. (55 i. No contribution 
proportional to the gravitational anomaly coefficient is found in this case. 



3.3 Conductivities for the Energy Flux 



We will include for completeness the result of the chiral magnetic and vortical con- 



ductivities for the energy flux, corresponding to the last two expressions in eq. (47 1. 

The chiral magnetic conductivity for energy flux, dj, follows from the correla- 
tion function of the energy momentum tensor and the current, and so it computes 



in the same way as the vortical conductivity in Sec. 3.1 From an evaluation of the 
corresponding Feynman diagram one finds that the result is the same as eq. ( [SO] ). 
Then one concludes that 

((j|)„ = ((Jv)„, (76) 



where ((Jv)a is given by eq. (73 i. Although these coefficients are equal, they de- 
scribe different transport phenomena. Whereas ((7j)u describes the generation of an 
energy flux due to an external magnetic field B„, {Gv)a describes the generation of 
the current due to an external field that sources the energy-momentum tensor r"'. 

Finally the chiral vortical conductivity for the energy flux, Gy, follows from the 
correlation function of two energy momentum tensors. There are three contributions 
out of the four possible terms. One of these terms involves a sum over fermionic 
Matsubara frequencies of the form 

^Y.{mf)^A,{mf ,q)A^{ia)f + i(0„, q + k) ^ ^{i(0„, Eg, Ey+k,t,u)+ (77) 
P a>f 



+ 



1 



ico„ + tEa - uE„+k 



tE'^n{Eg - tfl f) - u{Eg+k - uico„)^n{Ey+k - ujj.-^ 



where ^ corresponds to the ultraviolet divergent vacuum contribution which we 
choose to remove. The zero frequency, zero momentum, limit of the chiral vortical 
conductivity for the energy flux writes 



1 



a,b,c a 



(78) 



This coefficient describes the generation of an energy flux due to a vortex (or a 
gravito-magnetic field). The correlators (76 1 and (78 i enter the chiral magnetic and 
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vortical conductivities in the Landau frame, respectively, as defined in 112011211 122]| . 



see eqs. (49i-(50i. We have also checked that to lowest order in (O and k one has 
(j-OzjOz'^ _ where p is the pressure of a free gas of massless fermions, and 

(jOzyz) =0||26ll. 



3.4 Summary and specialization to the group i7( 1)^ x [/( 1)^ 

The results for the different conductivities are neatly summarized as 

= ^dabcil' , (79) 

1 

(C7v)« = (ctJ)„ = -T-^dabcll'' + l^ba , (80) 

(Ty - -^d,,b,ll"ll''ll' + '^ball" . (81) 
llK^ 12 

The axial and mixed gauge-gravitational anomaly coefficients are defined by 

date - \ [tr {Ta{T,„ r,})z, - tr (r„{7i, r,})fi] , (82) 

ba^^{Ta)L-^r{Ta)R, (83) 

where the subscripts L^R stand for the contributions of left-handed and right-handed 
fermions. The result shows that these conductivities are non-zero if and only if the 
theory features anomalies. 

For phenomenological reasons it is interesting to specialize these results to the 
symmetry group U{\)v x C/(1)a, i-C- one vector and one axial current with chemical 
potentials /i^ = ^ + ^1a, /Ir = ^-/Za, charges g'^^ = (1,1) and = (1,-1) for 
one left-handed and one right-handed fermion. We find (for a vector magnetic field) 

{(^bWv = ^ , {aB)Av = ^ , (84) 

iay)v = {al)v = ^ , iav)A = (oDa = + ^ , (85) 

=^(3M^ + Mi) + fr^ (86) 

Here {Ob)vv is the chiral magnetic conductivity ll25l . (Ob)/^ describes the gener- 
ation of an axial current due to a vector magnetic field [45], (c7v)y is the vector 
vortical conductivity, [Oy^A is the axial vortical conductivity, and Oy is the vortical 
conductivity for the energy flux. The vector and axial magnetic conductivities for 
energy flux ((Tj)y and {o^a coincide with the chiral vortical conductivities. 
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In this section for simplicity we will consider a holographic system which realize a 
single chiral symmetry with a gauge and mixed gauge-gravitational anomaly 
||28l . As we saw in the previous section in a more realistic model U{l)v x U{1)a 
the transport coefficients receive contribution from the gravitational part only in the 
axial sector. For a study of such a system with a pure gauge anomaly using Kubo 
formulae, see [35J . 



4.1 Notation and Holographic Anomalies 

Let us fix some conventions we will use in the Gravity Theory. We choose the five 
dimensional metric to be of signature (— , +, +, +, +). Five dimensional indices are 
denoted with upper case latin letters. The epsilon tensor has to be distinguished 
from the epsilon symbol by Eabcde = y/—g£{ABCDE). The symbol is defined by 
E{rtxyz) = +1. We assume the metric can be decomposed in ADM like way and 
define an outward pointing normal vector to the holographic boundary of an asymp- 
totically AdS space ha °^ g^^ with unit norm n^in'* = 1 . So that the induced metric 
takes the form 

hAB = gAB -riAriB- (87) 

In general a foliation of the space-time M with timelike surfaces defined through 
r{x) — const can be written as 

ds^ = {N^ +NAN^)dr^ + INAdx^dr + hABdx^dx'^ . (88) 

The Christoffel symbols, Riemann tensor and extrinsic curvature are given by 

^NP = \g'^^ [dNgKP + dpgKM - dKgNp) , (89) 

R'^ npq = dpr^Q - dgFj^p + r^r^Q - rgj^r^p , (90) 

c 1 

Kav = hA^cnv = ^£nhAB , (91) 

where £„ denotes the Lie derivative in direction of ha- Finally we can define our 
model. The action is given by 



1 



l6nG 



S = -^—r I d^xJ- 



M 



gMNPe^^^ (^ — pj^lppQjfJ^XR^BNpR^AQR 



-Sgh + Scsk. (92) 
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Sgh = ^ I d*xV-hK, (93) 

SKG JdM 

ScsK^~:^[ d^x^iXnMe'^™AtiKpLDQK'ii, (94) 
InG JdM 

where Sgh is the usual Gibbons-Hawking boundary term and Da is the induced 
CO variant derivative on the four dimensional surface. The second boundary term 
ScsK is introduced to reproduce the gravitational anomaly at general hypersurface. 

Lets study now the gauge symmetries of our model. We note that the action is 
diffeomorphism invariant, but they do depend however explicitly on the gauge con- 
nection Am. Under gauge transformations 5Am = ^ m^ they are therefore invariant 
only up to a boundary term. We have 

55 = /^^^/'x^/=/^^e'^'^^e« [^nMF,,pFQR + XnMR^ bnpR\qr] 

--^/ d^^hnMe^'^DN^Kpr^DQK'^. (95) 
AnG JdM 

Now without lost of generality we can choose the gauge N — I and A'a = 
which define the so called Gaussian normal coordinates, and the metric takes the 
form ds^ = dr^ + Yijdx'dx^ . After doing the decomposition in terms of surface in- 
duced and orthogonal fields, all the terms depending on the extrinsic curvature can- 
cel thanks to the contributions from Scsk^- The gauge variation of the action depends 
only on the intrinsic four dimensional curvature of the boundary and is given by 

8S = / d^x^e'""" i^F,„„Fu + Xr' j„„R^ ,„) . (96) 

l6%GJdM V3 'J 

This has to be interpreted as the anomalous variation of the effective quantum action 
of the dual field theory. As consequence of the discussion in the subsection 2.1 we 
can recognize the form of the consistent anomaly and use eq. (j9| to fix (f for a single 
fermion transforming under a t/(l)i symmetry. Similarly we can fix X by matching 
to the gravitational anomaly of a single left-handed fermion eq. ([22| and find 



K- _ 1 ^ _ ^ 

487rG ^ 967r2 ' ~ 167rG ~ 16%n^ 



(97) 



The bulk equations of motion are 



Gmn ~ AgMN — :^FmlFn^— ^F^gMN + '2XeipQg(^!^VB {f''^R^ n)^'^) , (98) 

V^F^^ = -e''^'Q''{KF^pFQp + XR^BNpR'AQR) ■ (99) 

A remarkable fact is that the mixed Chem-Simons term does not introduce new 
singularities into the on-shell action for any asymptotically At/S' solution, i.e. no new 
counterterm is needed to renormalize the theory. See II28J for a detailed discussion 
of the renormalization of the model and Appendix 1 to see the counterterms. 
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4.2 Applying Kubo Formulae and Linear Response 



In order to compute the conductivities under study using the Kubo formulae eq. ( 47 1, 
we will use tools of linear response theory. To do so we introduce metric and gauge 
fluctuations over a charged black hole background and use the AdS/CFT dictionary 
to compute the retarded propagators 1461 (471 . Therefore we split the backgrounds 
and fluctuations as, 

gMN = SmN + £ ^MN , (100) 

AM^A'i^^+eaM. (101) 

After the insertion of these fluctuations and background fields in the action and 
expanding up to second order in e we can read the on-shell boundary second order 
action which is needed to get the desired propagators 148 J . 



(102) 



where prime means derivative with respect to the radial coordinate, 0l is a vector 
constructed with the Fourier transformed components of and /imWj 

<t>'ir,xn = J (103) 

and jz/ and ^ are two matrices extracted from the boundary action and that we will 
show below. 

For a coupled system the holographic computation of the correlators consists in 
finding a maximal set of linearly independent solutions that satisfy infalling bound- 
ary conditions on the horizon and that source a single operator at the AdS boundary 
B6ll47ll48ll49 |. To do so we can construct a matrix of solutions F' j{k,r) such that 
each of its columns corresponds to one of the independent solutions and normalize 
it to the unit matrix at the boundary. Therefore, given a set of boundary values for 
the perturbations, (pj,, the bulk solutions are 

<Pl{r)=F'j{k,r)(p(. (104) 

Finally using this decomposition we obtain the matrix of retarded Green's functions 

G,j{k) = -21im {^,MiF''j{k,r)y + mj) ■ (105) 



The system of equations ( 98 1-( 99 1 admit the following exact background AdS 
Reissner-Nordstrom black-brane solution 
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As" = jj{-f(r)At^ + Ax^) + -^Ar\ (106) 

= 0(r)df = (^V-^^df, (107) 

where the horizon of the black hole is located at r = th and the blackening factor of 
the metric is 

The parameters M and Q of the RN black hole are related to the chemical poten- 
tial jX and the horizon rn by|^ 

M.| + | , e=^. (109) 

The Hawking temperature is given in terms of these black hole parameters as 

'-AiL?f^^-^- 2.ri ■ ''''' 

The pressure of the gauge theory is P = jg^^^a and its energy density is e = 3P due 
to the underlying conformal symmetry. 

To study the effect of anomalies we just turned on the shear sector (transverse 
momentum fluctuations) Oa and hf and set without loss of generality the momen- 
tum k in the y-direction at zero frequency, so a = x,z- Since we are interested in a 
hydrodynamical regime {k, (O ^ T), it is just necessary to find solutions up to first 
order in momentum. So that we expand the fields in terms of the dimensionless 
momentum p = k/AnT such as 

/i«(r)=/z|'""(r)+;,/z(')"(r), (111) 

with the gauge field redefined as Ba ~ Oa/ For convenience we redefine new 
parameters and radial coordinate 

A = — 2~ ; ^= — 2 — ' " = ^r^ ; (113) 

In this new radial coordinate the horizon sits at m = 1 and the AdS boundary at m = 0. 
At zero frequency the system of differential equations consist on four second order 



* The chemical potential is introduced as the energy needed to introduce an unit of charge from the 
boundary to behind the horizon A(oo) —A{rn) which corresponds to the prescription (B) in table^ 
Observe that we have left the source value A(<») = v as an arbitrary constant for reasons we will 
explain later. 
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equationsj^ The relevant physical boundary conditions on fields are: h"{0) = H", 
Ba{0) = Ba', where the 'tilde' parameters are the sources of the boundary operators. 
The second condition compatible with the ingoing one at the horizon is regularity 
for the gauge field and vanishing for the metric fluctuation |26 1. 

After solving the system perturbatively (see li28l for solutions), we can go back 
to the formula ( |105| l and compute the corresponding holographic Green's functions. 
If we consider the vector of fields to be 



(114) 



<pj (u) — {Bx{u),h^,{u),B^(u),h](u 
the and ^ matrices for that setup take the following form 



Diag 



-3a/, -3a/, - 
u u 



(115) 



^AdS+d 






-AKikii-<l>L^ 





-3a 

_ j_ 












\ 


—3a 



(116) 





/o 








\ 







3 



































(117) 



where = S^Ads+d + ^ct- [^Notice that there is no contribution to the matrices 
coming from the Chern-Simons gravity part, because the corresponding contribu- 
tions vanish at the boundary. These matrices and the perturbative solutions are the 
ingredients to compute the matrix of propagators. Undoing the vector field redefi- 
nition introduced in ( |1 \2\ the non-vanishing retarded correlation functions at zero 
frequency are then 



Gx.tx — Gz,tz. — 



Gx7 — —Gvx — 



V3Q 
4nGL^ ' 
iV3kQK ikvK 



InGrl^ 



6nG 



(118) 
(119) 



' The complete system of equations depending on frequency and momentum is showed in Ap- 
pendix 2. The system consists of six dynamical equations and two constraints. 
^ ^cr is coming from the counterterms of the theory. 
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G.xjz. — Gtx.z — ^G.jx — —Gtz.x — - — ^ A H -p, , (120) 

G,.,. = G,,,,= ^^^, (121) 
_ _ iV3kQ^ 4niV3kQT^ 

^tx,tz - -^tz,tx - H 6 1 TT^ • (i^^) 

2nGr^ Gr^ 



We can do an important remark observing eq. ( 119|l. Remember that we left the 



boundary value of the background gauge field eq. ( 107 1 arbitrary as a constant v. 
But as the U{1) symmetry is anomalous in the Field Theory side, physical quantities 
have to be sensitive to the source Aq, [^indeed as we can check they are. In particular 
if we choose the value v = which corresponds to formalism (A) in table [T] we 
need to include the counterterm eq. ([30]l in order to get the same propagator as at 
weak coupling. In fact in 13411351 it has been shown that in the case of a propagator 
between two vector currents, choosing this specific value for v the propagator would 
be zero, giving us in consequence a zero value for the chiral magnetic conductivity. 
Hence to be consistent with the scheme we are working at, let us just consider v 
as a source in the field theory. Therefore the real propagator is the one with v = 
because as is well known we have to set all sources to zero after taking the second 



functional derivative of the efective action. Finally using the Kubo formulae (47 1 we 
recover the vortical and axial-magnetic conductivities 

"^^""2^ Gr2 -l2^ + ~[2- ^^^^^ 

All these expressions coincide with the results in section [3] ( [79] l, ( [80| and (jSTJ if 
we specialize to dabc = 1 and ba = I- They are in perfect agreement with the lit- 
erature 1201 ED EH |26] except for the contribution coming from the gravitational 
anomaly which is manifest by the presence of the extra XT^. All the numerical co- 
efficients coincide precisely with the ones obtained at weak coupling; this we take 
as a strong hint that the anomalous conductivities are indeed completely determined 
by the anomalies and are not renormalized beyond one loop. Evidence for non- 
renormalization comes also from ifSOl where a holographic renormalization group 
running of the conductivities showed the same result at any value of the holographic 
cut-off. We also point out that the term that appears as undetermined integration 
constant in the hydrodynamic considerations in Il5l1 should make its appearance in 



' In principle Ao could be gauged away for the symmetric case and in consequence observables 
should not depend on its value. For example look at 1 35 1 to see how in presence ofaf/(l)vxJ/(l)/i 
symmetry with only the U{l)v conserved, propagators do not depend on the specitic value of the 
zero component of the vector gauge source Vb . 
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(Jy. We do not find any such term which is consistent with the argument that this 
term is absent due to CPT invariance. 

It is also interesting to write down the magnetic and vortical conductivites using 
eqs. ( [49] i and (BOb as they appear in the Landau frame to compare with the Son and 
Surowka form II22I 

~ SttGMLVZ + GM - 4k2 2 e+P J ^ ^^^'^> 

S, ^ 36^^ 2K?.THrl^-2L^Q^) _ ^i^ 2 n^i \ . 2n^i\ .,^7. 

~ AkGML^ GML^r^ ~ &n2 \^ 3 e+P I + 24 e+P J ' ^^^'> 

These expressions agree with the literature except for the XT^ term. A last comment 
can be done, the shear viscosity entropy ratio is not modified by the presence of the 
gravitational anomaly. We know that 77 0= limm^o ^ < T^'T^^' >k=Q, so we should 
solve the system at = for the fluctuations but the anomalous coefficients 
always appear with a momentum k (see Appendix 2), therefore if we switch off 
the momentum, the system looks precisely as the theory without anomalies. In [52 1 
it has been shown that the black hole entropy doesn't depend on the extra Chem- 
Simons term. 



5 Conclusion and Outlook 

In the presence of external sources for the energy momentum tensor and the currents, 
the anomaly is responsible for a non conservation of the latter. This is conveniently 
expressed through 142 J 



^^j, ^ ^,vp. (^I^^-^VA + ^^^Mv/^'^ ap.) , (128) 
where the axial and mixed gauge-gravitational anomaly coefficients, dabc and ba, 



are given by ( 23 1 and ( 24 1 respectively. 

We have discussed in Sec. |2] the constitutive relations and derived the Kubo for- 
mulae that allow the calculation of transport coefficients at first order in the hydro- 
dynamic expansion. We explained also subtleties in the definition of the chemical 
potential in the presence of anomalies. The explicit evaluation of these Kubo for- 
mulae in quantum field theory has been performed in Sec.[3]for the chiral magnetic, 
chiral vortical and energy flux conductivities of a relativistic fluid at weak coupling, 
and we found contributions proportional to the anomaly coefficients dn^c and ba- 
Non-zero values of these coefficients are a necessary and sufficient condition for the 
presence of anomalies ll42l . Therefore the non-vanishing values of the transport co- 
efficients have to be attributed to the presence of chiral and gravitational anomalies. 



For a four dimensional holographic model with gravitational Chern-Simons term and a scalar 
field this has also been shown in 1531 . 
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In order to perform the analysis at strong coupling via AdS/CFT methods, we 
have defined in Sec. |4]a holographic model implementing both type of anomalies 
via gauge and mixed gauge-gravitational Chern-Simons terms. We have computed 
the anomalous magnetic and vortical conductivities from a charged black hole back- 
ground and have found a non-vanishing vortical conductivity proportional to ^ T^. 
These terms are characteristic for the contribution of the gravitational anomaly and 
they even appear in an uncharged fluid. The behavior had appeared already pre- 
viously in neutrino physics ifTSl . In ll23l similar terms in the vortical conductivities 
have been argued for, but just in terms of undetermined integration constants with- 
out any relation to the gravitational anomaly. Very recently a generalization of the 



results (123l-(125i to any even space-time dimension as a polynomial in ji and 
T [29 1 has been proposed. Finally, the consequences of this anomaly in hydrody- 
namics have been studied using a group theoretic approach, which seems to suggest 
that their effects could be present even atT—Q |54|. The numerical values of the 
anomalous conductivities at strong coupling are in perfect agreement with weak 
coupling calculations, and this suggests the existence of a non-renormalization the- 
orem including the contributions from the gravitational anomaly. 

There are important phenomenological consequencies of the present study to 
heavy ion physics. In ||55]| enhanced production of high spin hadrons (especially 
Q~ baryons) perpendicular to the reaction plane in heavy ion collisions has been 
proposed as an observational signature for the chiral separation effect. Three sources 
of chiral separation have been identified: the anomaly in vacuum, the magnetic and 
the vortical conductivities of the axial current j'^. Of these the contribution of the 
vortical effect was judged to be subleading by a relative factor of 10 . The T term 



124| i leads however to a significant enhancement. If we take ji to be the baryon 



chemical potential /i sa 10 MeV, neglect jx^ as in |55j and take a typical RHIC 
temperature ofT — 350 MeV, we see that the temperature enhances the axial chiral 
vortical conductivity by a factor of the order of 10^. We expect the enhancement at 
the LHC to be even higher due to the higher temperature. 

In this review we have presented the computation of the transport coefficients, 
and in particular their gravitational anomaly contributions, via Kubo formulae. It 
would be interesting to calculate directly the constitutive relations of the hydrody- 
namics of anomalous currents via the fluid/gravity correspondence within the holo- 
graphic model of Sec. [4] ||20ll2ni56l . This approach will allow us to compute trans- 
port coefficients at higher orders Il57ll58l . This study is currently in progress | 
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The result one gets for the counterterm coming from the regularization of the bound- 
ary action of the holographic model in section]?] is 



Sc, 



^TlG JdM 



1 1 

^ + 2^-12 



p'pJ - 
J 1 



where hat on the fields means the induced field on the cut-off surface and 



(129) 



(130) 



As a remarkable fact there is no contribution in the counterterm coming from the 
gauge-gravitational Chern-Simons term. This has also been derived in ll60l in a sim- 
ilar model that does however not contain Scsk- 



Appendix 2: Equations of motion for the shear sector 



These are the complete linearized set of six dynamical equations of motion, 



/(«) 
3 



= fi'a (m) + '^Ba («) + -if^ {cO^-f{u)k^) Ba [u) - 

UJ [Uj 



hf{u) 



\ uf{u) V3a / f{u 

= hf"(u) - ''^^ - (k^h"iu)+h" (u) (Ok) - 3auB'(u) 

u uf\u) 



(131) 



(2W -6(1 -f{u))) + (W -6(1 - f{u)))B'p{u) 



+2u{uhf' (u))' (ji^,{u)(0k + hf(u)k^ 



ua"r..\ , {[M^ua't 
f/u 



= (m) + (m) + -^^^^2 i^^'^fi'^) + 0}khf{u)) +2uikle, 

hi'iu) b^ 



+ (9/(m) - 6 + 3au') + ^ (^cokhf (m) + w^h^ (u) 



/(«) /(«) 
and two constraints for the fluctuations at (O,ky^0 



(132) 
uh^ (u) 
(133) 



= 0) [h" (m) - 3auBa{u)j +f{u)kh" (u) + ikXtap 
+ {9au^ -6{l-f{u)))Bp{u)] . 



lu^ {ahf +f{u)kh^.' {u) 

(134) 
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